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A MODIFIED MULTHOPP APPROACH FOR PREDICTING 
LIFTING PRESSURES AND CAMBER SHAPE FOR 
COMPOSITE PLANFORMS IN SUBSONIC FLOW 

By John E. Lamar 
Langley Research Center 

SUMMARY 

This report presents a modified version of Multhopp's subsonic lifting- surface 
theory which has been programed in two parts for the IBM 7094 electronic data processing 
system or Control Data 6400 computer system along with a discussion of the character of 
its results. The first part is used to find both basic and additional loadings over a given 
planform with known mean camber surface and the second part is used to determine the 
required mean camber surface for a given planform and set of loadings. 

For the loading program, various aerodynamic characteristics are determined on 
both simple and composite planforms when the spanwise loading is symmetrical. Studies 
are conducted to determine when these answers are most valid, and some results for 
delta, sweptback and tapered, double delta, and variable-sweep wings are compared with 
other theories and experiments to determine the accuracy of this method. 

Application of this method is then made in predicting the aerodynamic effects of 
changing the outer-panel sweep of a variable- sweep wing, of increasing the Mach number 
in the subsonic regime, and of incorporating twist and camber in a wing. The first two 
applications also have experimental data presented for comparison. 

A study is also undertaken for the mean camber surface program to determine when 
its best results are obtained. Subsequently, comparisons between the present and other 
methods are made for a two- and three-dimensional case to aid in the evaluation of the 
present method. One application is made for a highly sweptback and tapered planform. 

INTRODUCTION 

In the process of wing design, one is required to determine either the load distribu- 
tion from a given planform and mean camber surface or the mean camber surface from a 
prescribed planform and loading. From previous investigations a considerable amount of 
information of both a theoretical and experimental nature is available to aid the designer 
if the wings have delta or simple sweptback planforms. However, planforms now being 



considered (which in many cases have the divergent requirements of supersonic or hyper- 
sonic cruise and subsonic loiter or ferry capability) are much more complex, involving, 
in many instances, cranked leading and trailing edges. Planforms of this type are classi- 
fied as being of a composite arrangement. 

For these planforms not as much general information is available, due both to their 
newness and to the large number of planform variables. First-order solutions to these 
problems have been found in recent years for the sonic (refs. 1 and 2) and supersonic 
(refs. 3 and 4) speed regimes which are applicable to arbitrary wing planforms. Refer- 
ences 3 and 4 have been programed for the high-speed electronic computer so that the 
designer does not have to perform the lengthy mathematical manipulations required to 
obtain a solution for each wing considered. 

In the subsonic speed regime, the basic theoretical method of Multhopp (ref. 5) has 
received wide acceptance as one of the most accurate methods for predicting aerodynamic 
loading data; consequently, it has been selected as the approach to use in solving these 
two design problems for composite planforms. Some reasons for this selection are 

(1) The chordwise pressure distributions (and other section data) as well as overall 
aerodynamic characteristics are determined 

(2) The spanwise locations of the lift-producing singularities and control points are 
more concentrated along chordwise rows near the tip to insure an adequate representation 
of the spanwise load distribution 

(3) The method can as Multhopp mentioned in reference 5 be used for wings with 
planform kinks at locations other than just at the plane of symmetry 

A computerized solution of the loading determination problem for arbitrary wings 
has been developed from the basic Multhopp method by Van Spiegel and Wouters in refer- 
ence 6 and is referred to as the modified Multhopp method. Certain changes in and exten- 
sions to their work have been made in the present report to increase the accuracy in 
representing the actual chordwise pressure distribution and finding the influence of each 
of these distributions on each point where the boundary conditions are to be met. These 
modifications resulted from adding additional pressure modes (up to 10 may be used) and 
by replacing the polynomial approximation technique used in the chordwise integration 
with the Gaussian quadrature method. 

In addition, the formulation procedure and all the modifications made to the basic 
Multhopp approach are described as well as the techniques used in its solution for both 
the load and mean camber surface determination problems. After this a study of the 
resulting answers is made to determine their sensitivity to the number and location of 
the points where the boundary conditions are met or calculated. Then applications of the 
programs are made to an assortment of wings, including composite wings and the results 
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are compared with other theories and, where possible, experiment. All the predictions 
made from this modified Multhopp method have been obtained with the aid of the IBM 7094 
electronic data processing system and Control Data 6000 series electronic computers. 

A listing of the computer programs used to find the surface loadings (Langley pro- 
gram A0313) and to find the mean camber surface (Langley program A0457) is pre- 
sented along with a description of the input data required, sample input and output data, 
and pertinent comments in a "Supplement to NASA TN D-4427." A large portion of these 
two programs is given over to the computation of the geometric representation required 
by later parts of the program from the input quantities supplied. Two additional programs 
(Langley programs A1590 and A1591) which are useful in obtaining geometric input data 
for either program A0313 or A0457 are also presented in the supplement along with brief 
descriptions and sample cases of each. A request form is included at the back of this 
paper. 


SYMBOLS 


Any convenient system of measure can be used as long as the linear dimensions of 
the planform are based on a semispan of unity. 

A aspect ratio, b2/S 


a.c. 


b 


aerodynamic center, in fractions of 


erence chord (positive aft), 


-ec 


m 


9C L 


c ref > referenced to leading edge of ref- 


+ 


1 

4 


wing span, set equal to 2 




Multhopp’ s integrating coefficients 


Ca 

Cb,o 

c D,i 

c D,ii 

c L 


axial-force coefficient, Axial force 

*loo®ref 

b f 1 

root bending- moment coefficient, = \ ( c ^ c )basic^ d? 7 

I* 61 Q 

induced drag based on the spanwise distribution of circulation 

induced drag based on axial force, aC l - 

lift coefficient, Wmg lift 
9ooSref 
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overall lift -curve slope 


C L, 

C 


'a 


m 


'm 


a 


ACr 


Cs 


. , Wing pitching moment c ref 

pitching-moment coefficient, about — - — 

QocSrefCref 4 

overall pitching-moment- < 


Cpgf 

-curve slope about — — 


incremental pressure coefficient, 


Pupper ~ Plower _ Ap 

*00 " q - 


leading- edge -suction coefficient for zero degrees of leading- edge sweep, 


fab r 


r 1 1 


J_1 c(?7) 

qoo 


dri 


streamwise chord at y ,77 


mean geometric chord, 


f b/ f o2dy 

J -b/2 

r b/2 

\ c dy 
J -b/2 


c(i?),c(y) streamwise half-chord at n,y 
average chord, S/b 


'av 


c l 

c m 


a 


c n, 


a 


c ref 


local lift coefficient, 


Local chord load 
qccC 


local pitching- moment- curve slope about local leading edge 
local slope of normal-force coefficient 
root chord 

reference chord, may be mean geometric chord of total wing c except 
when the planform has an inboard trailing- edge chord- extension or may be 
mean geometric chord of only wing outer panel extended to plane of symme- 
try in streamwise tip position 


d(r?) 


x- location of mid chord at 77 
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f j 

K(4>,y;3,Tj) 

Kj(^,y;r?) 

Lyn^s) 

M 

m 

N 

Ap 

qtjfa) 

Qoo 

s 

Sref 

s 

U 

w 

x,y,z 


jth chordal loading function 
subsonic kernel function 

subsonic kernel function integrated with jth loading across the chord 

influence of jth loading function at r/ n on control point at y v and 
x = (c/2)(l - cos (p s ) 

elements of influence matrix when n 4- v 

Mach number 

number of span stations where the pressure modes are defined 
number of chordal control points at each of m span stations 
lifting pressure, p upper - p lower 

coefficient of jth chordal loading function for additional load 
coefficient of jth chordal loading function for twist and/or camber load 
dynamic pressure 
total wing area 

reference area, may be either total wing area or wing area of wing outer 
panels extended to plane of symmetry in streamwise tip position 

used to locate chordal control points and ranges from 1 to N 

free- stream velocity 

perturbation velocity in z-direction 

rectangular Cartesian coordinates nondimensionalized with respect to b/2 
where origin is in plane of symmetry at half root chord (fig. 1) (they are 
associated with control points) 


5 


x c 


x c r ef 


^cp’ x cp 
a 


®i 


a l 


x-location of leading edge of c 

x- location of leading edge of c re f 

centers of pressure 

angle of attack, deg 

induced angle of attack, radians 

local angle of attack, radians 

& ’ dx U 


j3 Prandtl-Glauert compressibility factor, \Jl - 

y local nondimensional lift coefficient, c i c /Zb 

A wing tip skew angle, deg 

Q angle for locating control points along span (see fig. 2(a)) 

d angle used to locate pressure doublets chordwise, 0 at leading edge and fl- 

at trailing edge 

k ratio of distance of leading- edge break from plane of symmetry to b/2 

A outboard leading- edge sweep angle, deg 

Tip chord 

X taper ratio, root chord 

H ratio of distance of trailing- edge break from plane of symmetry to b/2 

rectangular Cartesian coordinates nondimensionalized with respect to b/2 
along X- and Y-axes, respectively, and associated with pressure doublet 
locations 

<p angle used to locate chordwise rows of control points (see fig. 2(b)) 

X inboard leading- edge sweep angle, deg 
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inboard trailing- edge sweep angle, deg 


H' 

Cl outboard trailing- edge sweep angle, deg 


Subscripts : 

cp 

d 

le 

n,v 

0 

P 

te 

x 

y 

1 


center of pressure 
distributed 
leading edge 

suffixes numerating spanwise stations, v station being influenced and n 

station doing the influencing; n ,v range from - m ~ . . . 0, . . . 21_ 

2 ^ 

value taken at Cl = 0 
pivot 

trailing edge 

x-direction 

y-direction 

result due to twist and/or cambered wing at zero root chord angle of attack 
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BASIC FORMULATION 


The determination of surface loadings on both plane and warped wings and the deter- 
mination of the mean camber surface can be made from a number of theories already in 
use. All the theories share a common feature in that the wing is replaced by a mathemat- 
ical model which produces a similar downwash field to that of the wing. However, the 
methods differ mainly in the potential used to produce the representative perturbations of 
the flow field around the wing. These potentials are either of the velocity or acceleration 
type (ref. 7). The original Multhopp method and its present modified version use the 
acceleration potential approach in conjunction with the linearized Euler equations to relate 
the pressure difference across the wing to the downwash field on the wing surface. This 
results in the following equation which is derived in detail in reference 6 for the incom- 
pressible case: 
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( 1 ) 



i r 1 i 

8 ^ J-l*. 


rke Ap (|, t?) 

J «je ( y -^) 2 3 4 


1 + 


(X - |)//3 


/(x - ?) 2 


Ml d?7 


(y - ^) 2 


Reduction of Basic Equation to Matrix Form 

Equation (1) is difficult to solve for the local pressures Ap(|,?j) since they appear 
as part of the integrand and because of the integration across the second-order singular- 
ity. The term w/U and its resulting mean camber surface is also difficult to deter- 
mine because of this same required integration. A solution to either problem can, how- 
ever, be made by 

(1) Replacing Ap (4,r?) with as many of the 10 pressure modes of unknown or known 
amplitude as are needed 


Ap(0,i?) 


%(v) 

c(t?) 


cot 


i/ 

1*1 

j=i 


c(v) 


sin jO 


( 2 ) 


where 

| = -c(t))cos d + d(rj) (3) 

(2) Employing the Multhopp quadrature formula to approximate the spanwise inte- 
gration of the pressure doublets; this approximating procedure has another feature in that 
every spanwise station need not be represented but only those stations where information 
is required 

(3) Using a logarithmic singularity correction term (appendix A) to account for the 
integration across the second-order singularity (this amounts to finding the finite part) 

(4) Setting up equation (1) with these approximations in a matrix equation of the form 



The details of formulating the matrix of influence coefficients are given in appen- 
dix A where they have been taken for the most part directly from reference 6 and pre- 
sented here for the reader's convenience. 


The points (determined from 0 S and y v ) which are selected either to specify the 
tangential flow boundary condition or to find the downwash ratio are called control points 
and are located in both the spanwise and chordwise directions by the use of semicircles 
which are laid along the span from tip to tip and along the local chord from leading to 
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trailing edge. The spanwise positions, which result from an equiarc angular displace- 
ment about the plane of symmetry, can be represented by the following equation: 


Yv ~ sin 


m + 1 


( 5 ) 


where m is the total number of spanwise stations from tip to tip and v ranges from 
- m 2 ~ 1 to m 2 1 in increments of 1. (See fig. 2(a).) These positions are those 
required by the Multhopp quadrature formula for the integrand evaluation when approxi- 
mating the value of the spanwise integral (ref. 5). 


At each spanwise location, the control points are also located for equiarc angular 
displacements defined by 


0s - 


27TS 

2N + 1 


( 6 ) 


where N is the total number of chordwise control points and ranges from 1 to N. (See 
fig. 2(b).) Note that the control points never reach the leading or trailing edge. From 
this, the control-point locations are found to be determinable by the following equation: 


| = |(1 - cos 0 S ) 

which locates them aft of the local leading edge. (See fig. 3.) 


(7) 


In order to solve for either the loading (inverse problem) or mean camber surface 
(direct problem), equation (4) is used. The loadings are found by (1) inverting the matrix 

qj fan) 


in equation (4) and thus solving for the 


array and then (2) using these values in 


the equations given in appendix B to determine the local and overall aerodynamic loading 
description. The mean camber surface is found by the direct matrix multiplication of 
equation (4) along with the procedure described in the mean camber surface determination 
section of this report. 


Modification Made to the Original Multhopp Method 

From a study of Multhopp’ s original method, several features are seen to need 
modification in order that a computerized solution of it could be made to the best advan- 
tage. Aside from the logarithmic singularity correction term whose best form is dis- 
cussed in appendix A they are as follows : 

1. The method for solving for the coefficients of the loading terms in the inverse 
problem. In Multhopp’s report, a scheme of iteration on the downwash is used, this is 
improved in reference 8 with a scheme of submatrix inversion and in reference 6 with a 
scheme of whole matrix inversion which is also used herein. 
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2. The chordal integration that sums the influence of the pressure doublets on con- 
trol points not at the same span station. In Multhopp’s report, it is developed from a 
MacLaurin series then graphed for use. Van Spiegel and Wouters (ref. 6) follow the 
scheme developed by Van de Vooren (ref. 9) with which the value of the integral is deter- 
mined, by (a) approximating the loading functions with polynomials, (b) evaluating them at 
discrete locations, and (c) multiplying the results of part (b) with predetermined coeffi- 
cients. In this report, the polynomial approximation of reference 6 is replaced with a 
Gaussian 20- to 30-point numerical integration scheme. 

3. The number of chordal loading functions and, therefore, control points. In 
Multhopp's report they are limited to two, cot which gives lift with no pitching 
moment about c/4 and cot k ~ 2 sin 3 which yields no lift but a pitching moment about 
c/4. From these, the known downwash and the matrix of influence coefficients, one is 
able to solve for the coefficients of the loading terms which are the local circulation and 
pitching moment about local c/4. Garner and Lehrian (ref. 10) have extended to four the 
number of chordal loadings (and consequently, control points) that may be used but employ 
combinations of terms rather than separate loading modes. Van Spiegel and Wouters also 
extend the number of chordal loadings and control points to four but separate the loadings 
into cot sin 3, sin 25, and sin 35, whose coefficients are the strength of the local 
pressure doublets supporting that loading contribution over the chord. In this report 
they have been extended to 10 with the addition of sin 45, sin 55, sin 65, sin 75, 

sin 85, and sin 95 for additional flexibility. 

Multhopp’s method with the aforementioned modifications is hereafter referred to 
as the present method and has been programed to handle a wide variety of composite 
planforms, with a break in the leading and trailing edges as well as the variable- sweep 
wings. The program at present is capable of handling sets of chordwise and total span- 
wise stations of 1 by 41, 2 by 41, 3 by 41, 4 by 41, 5 by 39, 6 by 31, 7 by 27, 8 by 23, 9 by 
21, or 10 by 19, with the number of chord stations given first. This is not a theoretical 
limitation but is near the upper limit of machine storage presently available. Also, with 
the sets of stations listed above, almost any planform of low to moderate aspect ratio can 
be handled. 

The distribution of lift obtained by the present method is dependent on the subsonic 
Mach number (see appendix A) at which it is computed. In this report, the Prandtl- 
Glauert relations (ref. 11) have been used to account for this change but are valid only as 
long as ACp varies linearly with a. 
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LOAD DETERMINATION 


Because of the nature of the numerical solution, different combinations of N and 
m for the same planform can result in different answers for the aerodynamic character- 
istics. Hence, a study was undertaken to try and optimize their combination so that rea- 
sonably accurate results are obtained. The results of this study are presented in the next 
section. Then comparisons of other theoretical methods and experiment with the present 
method are made to gain confidence in the answers obtained with it. And lastly, illustra- 
tive examples are given to show some applications for which the present method can be 
used. 


Optimization of N and m 

For any given planform and pair of N and m, a set of chordal loading coefficients 
qj (r] n ) are obtained from the solution of the inverse problem for the flat plate and are 
used to determine the loading characteristics of the planform. However, since a different 
set of qj(?7 n ), and therefore loadings, will be found for each pair of N and m, some 
criterion is needed by which the most appropriate set of qj(? 7 n ) can be selected. One 
criterion, suggested by Multhopp (ref. 5), is to examine the chordwise load distribution by 
paying particular attention to the leading- edge suction force developed. The contention 
being that when the overall leading- edge suction force, which is assumed equal to the 
axial force and given in coefficient form for an unswept planform by 


c s = 


j37rb r 


8S re f ^-1 c(t?) 


%(v) 


n 2 


dri 


is such as to cause the quantity 


( 8 ) 


c D,ii “ QiC L - C A 

to be equal to the induced drag coefficient given by the equation 


C D ,i 


b2 r 1 
Sref J - 1 


ya.i dr] 


( 9 ) 


( 10 ) 


then the cotangent loading (eq. (2)) can be considered valid. This contention is based on 
the conclusions of a study made by Munk in reference 12 where he found that the term 
c D,i was independent of the chordwise loading distribution. In the discussion that fol- 
lows, the result of equation (9) is called the suction value and the result of equation (10) 
is called the circulation value. 
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Upon comparing graphically the two coefficients just discussed for many combina- 
tions of N and m for two rectangular planforms of different aspect ratio, as is done 
in figures 4(a) and 5(a), it is noted that whereas the circulation value is essentially invari- 
ant with N and m, and is plotted as a straight line without symbols, the suction value 
shows a marked dependence upon N and m. If values for the aerodynamic characteris- 
tics at each pair of N and m where the two solutions are first equal or have a mini- 
mum difference are plotted as solid symbols, it is seen that the aerodynamic center and 
CLq, values for the A = 2 planform (figs. 4(b) and 4(c)) do not vary appreciably from 
results obtained when N = 8 at the higher values of m. However, in the aerodynamic 
center and Cl^ graphs (figs. 5(b) and 5(c)) for the A = 7 planform, it appears that the 
values predicted do not necessarily agree with those at the higher N and m. One rea- 
son suggested for the disagreement is that not enough values of m were available to be 
used at the higher N values (storage limitation) so that the resulting answers found by 
convergence could have been established. Note that the aerodynamic center is the more 
sensitive of the two; hence, a major part of the discussion that follows is concerned with 
it. 

From a study of figures 5(b) and 5(c), the curves for N = 2 and 4 seem to show 
that a particular relationship must exist between N and m in order for the answers to 
reach a level of convergence. Whatever the relationship is, it leads to results which are 
close to those given by the solid symbols. From the N = 4 curve in figure 5(b), it 
seems that a good result from the point of view of convergence would occur at m = 37, a 
solid symbol. Now if this relationship between m and N is valid, then its multiplica- 
tion with lower N values to determine new m numbers should lead to results which 
are as good as those found from the original N and m values. 

An application of this procedure would be that of multiplying the ratio of (37/4) by 
an N value of 2, thereby yielding an m number of approximately 19. When the curves 
of aerodynamic center and Cl^ at N = 2 and m = 19 are examined in figures 5(b) 
and 5(c), this combination of N and m is seen to give results which approximate well 
those obtained for N = 4 and m = 37. This occurs even though N = 2 and m = 19 is 
not a solid symbol which indicates that even though the ratio was obtained from a combi- 
nation of N and m that was a solid symbol, applications of it at lower values of N 
do not necessarily result in combinations of N and m which are other solid symbols. 

It should be noted that had the N = 2 and m = 13 solid symbol been selected as the 
numbers to be used in determining the ratio of m to N, the m value which would have 
resulted from the multiplication of the ratio with N = 4 would be 26. This combination 
of N and m does not give answers which approximate as well the results found from 
convergence. 
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Hence, sets of N and m can be calculated that will result in answers for aero- 
dynamic center and Cl^ which approximate well those found at a solid symbol provided 
they have been obtained from a solid symbol which has a larger N value than for the 
combination sought. In addition, from the results of both figures 4 and 5, the relationship 
between N and m is seen to indicate that an aspect ratio dependency also exists, 
requiring the use of more m for a given N at the higher aspect ratios. This is dis- 
cussed in more detail later. 

Examining again the variation of the aerodynamic center with N and m for the 
A = 2 planform in figure 4(b), it is noted that a part of the aerodynamic- center movement 
can be correlated with the suction force coefficient which is seen implicitly in figure 4(a). 
In general, as the suction increases ^a decrease in CD,ii/CL^) a corresponding forward 
movement in the aerodynamic center is observed. This simple relationship implies that 
if any two sets of N and m give rise to the same suction coefficient, then their aerody- 
namic centers will also be the same. However, this is not generally true. Look, for 
example, at a value of Cd,u/Cl^ of 0.20 and note all of the different combinations of 
N and m that will give it and then the corresponding aerodynamic centers. Some dif- 
ferences are seen to exist between the aerodynamic centers predicted by the different 
combinations of N and m. 

There are several possible reasons for this difference. The first one is that the 
relationship assumed between the leading- edge suction and the aerodynamic center is not 
a completely valid one since it depends not only on the qg(??n) co t ^ terms (which give 
rise to the leading-edge suction) but also the q 2 (? 7 n )sin 23 terms. 

The second reason is that since in the solution for the qj(r 7 n ) set the number of 
control points is required to be equal to the number of unknown coefficients, the flat plate 
representation of the planform may not be made at enough control points. This could be 
remedied by an overdetermined solution (i.e., one where there are more control points 
and number of equations than unknown coefficients). 

The third reason follows from the second in that, with the boundary condition only 
specified at a maximum of 10 points chordwise, the variation of w/U that results by 
matrix multiplication using the q j (^n.) set determined from the specified boundary con- 
ditions does not assume the boundary condition value except at the control points. Else- 
where it is either higher or lower and because of the logarithmic singularity correction 
term (eq. (A9)) it is singular at the leading and trailing edges. This variation in w/U 
when acted upon by the chord loading distribution should also give rise to an axial force 
heretofore unaccounted for. It is called the distributed axial force and is computed by 



Numerical results have not been obtained for c j since every point along the chord at 
each spanwise position must contain the influence of every loading function at every other 
spanwise location in order that the quantity w/U can be determined. This implies 
evaluating integrals similar to those used in finding qj(%) but at an extremely large 
number of locations. In addition, the singularities at the edges of the chord would be dif- 
ficult to evaluate without a coordinate transformation. (Note that if the mathematical pro- 
file were truly flat, then all the axial force generated would be concentrated at the leading 
edge in the form of a suction force, since w/U would be everywhere equal to a.) Thus, 
only one solution for the axial-force coefficient has been obtained and is derived solely 
from the leading- edge suction. 

From figures 4 and 5, it appears (as one would probably expect) that for the most 
part at the higher N and m values the answers for and aerodynamic center are 

in what appear to be converging situations. However, it has also been found that one need 
not always choose the maximum N and m to achieve these answers but may by proper 
insight and judicious selection obtain equally as valid answers at lower values of N and 
m. Thus for planforms which have zero sweep on the leading edge, it is suggested that 
pairs of N and m be selected for the particular planform at the desired Mach number 
in a manner so as to emphasize its largest dimension; that is, planforms having large 
values of j3A would require a pair of N and m where m was optimized after a 
suitable value for N was established, and planforms having low values of (3A would 
require that a large value for N be used with a corresponding reduction in m. 

Selecting pairs of N and m by this procedure should give answers that will be in or 
be very close to those in the region of converging results. This can be accomplished 
approximately by using the following expression: 

m = (4 to 5)j3A^ (12) 

A similar type of analysis which leads to equation (12) should be made for wings with 
sweepback using the leading-edge-suction equation recently developed by Garner in refer- 
ence 13. 

In an attempt to see if the suggestions given would be applicable to arbitrary plan- 
forms, aerodynamic-center and Cl graphs were prepared for two double delta plan- 
forms which are shown in figures 6 and 7. It appears that even with taking N = 8 the 
formula does not give quite enough m (only 20) to predict the aerodynamic center at the 
higher m values (figs. 6(a) and 7(a)). Consequently, for low-aspect-ratio wings with 
sweepback, an N value of from 6 to 10 should be used with as many m values as the 
machine storage limitation will allow. 

In addition to the overall aerodynamic characteristics which have already been dis- 
cussed, the influence of different combinations of N and m has been examined with 
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regard to their influence on the span load distribution. (See figs. 8 to 11.) From this 
examination, the reason becomes apparent why the induced drag parameter, 
was stated to be essentially invariant with different combinations of N and m. The 
reason is that the span loadings themselves appear for the most part to be independent of 
N and m. 


Comparisons and Illustrative Calculations 

The following section contains a number of comparisons with other theoretical 
methods and experimental data which have been used to verify the results obtained with 
the present method. After the verifications have been made, illustrative calculations are 
made to show the present methods' usefulness in computing sweep, Mach number, twist, 
and camber effects. Additional illustrative calculations of the present method have been 
made and are given in reference 14. 

Comparisons wit h other theoretical methods .- In order to insure that the computa- 
tional procedures presently employed give equivalent results to those obtained by Multhopp 
(ref. 5) and Van Spiegel and Wouters (ref. 6), table I has been prepared to give a tabular 
comparison of CLq,, aerodynamic center, and CD,i/fc;L 2 - The indications from table I 
are that all three methods give results for C;^ and CD,i/t;L^ within 2 percent of 
each other and the aerodynamic centers are within 0.3 percent of the reference chord of 
each other when at a constant value of N. When the values are compared 

with 1 /ttA (0.07457 for this case) they are found not to differ from it by more than 
2.2 percent, indicating that the computed span loading distributions are nearly elliptical. 

With the present method established as giving reliable results when compared with 
the basic and modified Multhopp procedures (refs. 5 and 6) and with the studies previously 
conducted, which indicated that an extension in the number of chordwise pressure modes 
was needed for certain planforms, comparisons of predicted results obtained with the 
present and other methods are made. The following planforms are offered as examples 
of this : 

(1) The aerodynamic centers and CLq, values for the four simple planforms shown 
in figure 12 are compared in tables II and m, respectively, with results obtained by 
DeYoung and Harper in reference 15 by using a Weissinger seven-point solution. (Values 
of m of 9 or 11 were chosen since N = 4, A = 2, and /3 = 1.) 

From table n the aerodynamic centers predicted by the present method are seen to 
lie from 1.7 to 4.2 percent of the reference chord ahead of those given by reference 15 
and from table m the values of CLq, are seen to vary from 1.62 to 5.35 percent below 
those computed by using the present method. These differences are attributed to the 
more accurate accounting of the induced camber effects by the present method. 
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TABLE I.- COMPARISONS OF RESULTS OBTAINED FOR THIS PLANFORM 
WITH THE PRESENT AND THE MULTHOPP METHODS 
OF REFERENCES 5 AND 6 



i 

.300 

* 


Method 

N 

m 

C 

a.c. 

c D,i/c L 2 

Ref. 5 

1 

15 

0.0564 

0.270 

0.08125 

Present 

1 

15 

.0569 

.271 

.08129 

Ref. 5 

2 

15 

0.0571 

0.279 

0.08069 

Ref. 6 

2 

15 

.0580 

.278 : 

.08093 

Present 

2 

15 

.0580 

.277 

.08087 

Ref. 6 

3 

15 

0.05836 

0.262 

0.08090 

Present 

3 

15 ( 

.05834 

.263 

.08086 


TABLE H.- AERODYNAMIC CENTERS FOR 
A = 2 PLANFORMS OF FIGURE 12(a) 
WHEN N = 4 



Aerodynamic center for - 

Method 

A = 0°, 
X = 1.0, 
m = 11 

A = 0°, 

X =0.25, 
m = 11 

A = 45°, 
X = 1.0, 
m = 9 

A =45°, 
X =0.25, 
m = 9 

Present 

0.208 

0.191 

0.156 

0.248 

Ref. 15 

.250 

.245 

.192 

.265 


TABLE m.- LIFT-CURVE SLOPE FOR 
A = 2 PLANFORMS OF FIGURE 12(b) 
WHEN N = 4 



Lift -curve slope for 

- 

Method 

A = 0° 

X = 1.0, 
m = 11 

A = 0°, 

X =0.25, 
m = 11 

A = 45°, 
X = 1.0, 
m = 9 

A = 45°, 
X =0.25, 
m = 9 

Present 
Ref. 15 

0.0432 

.0422 

0.0427 

.0415 

0.0398 

.0376 

0.0432 

.0425 
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TABLE IV.- DOUBLE-DELTA NONDIMENSIONAL LIFT COMPARISON 






C0S ( A c/2) eff 

A 

CL a 

C L« 

K 

Aspect ratio 

X 

cos ( A c/2) eff 

A 

(Present) 

(a) 

A 

(Ref. 16) 

0 

1.88810 

0 

0.7559 

2.4979 

0.02086 

0.02108 

.25 

1.41500 

82.63 

.5212 

2.7148 

.02234 

.02040 

.50 

1.24920 | 

78.07 

.4800 

2.6025 

.02253 

.02075 

.75 

1.11830 

75.00 

.5939 

1.8830 

.02309 

.02309 

1.00 

3.90020 

40.00 

.9221 

4.2297 

.01545 

.01641 


a N = 8; m = 23. 


Comparisons of the span loadings are also made between the present method and 
reference 15 for two of the planforms in figure 12 and are presented in figures 8 and 9 
for various combinations of N and m. The span loading distributions computed by the 
two methods are in good agreement for most sets of N and m. 

(2) A series of double-delta planforms, whose general planform parameters are 

given in table IV, have been used to obtain values of CLq^A which are also presented in 
table IV. They are compared with the values of Ci_, a /A predicted by Spencer in refer- 
ence 16. From table IV, the predicted values show that a maximum difference 

of 6 percent or -9 percent exists between the results of the two methods for the planforms 
examined. The difference is attributed to the method used in reference 16 because it does 
not consider the distribution in lifting pressure over the surface as the present method 
does but is a refinement of a mathematical expression relating A and sweep to the value 
of Cl for wings which have broken leading edges. 

(3) The three predicted span loading curves, given in figure 13 for a sweptback and 
tapered planform, come from references 15 and 17 and the present method. The span 
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loadings are seen to differ, in terms of y cp , within 0.9 percent of the semispan from each 
other. The Ct results of the present method and Campbell’s method (ref. 17) agree 

■ Lj G! 

quite well and are about 5 percent greater than the result of DeYoung and Harper (ref. 15). 
For the aerodynamic center, the results of the present method and Campbell's method 
(ref. 17) are 3.2 and 2.1 percent, respectively, of the reference chord ahead of that pre- 
sented by DeYoung and Harper (ref. 15). 

(4) The two span loading curves (obtained from ref. 8 and the present method) shown 
in figure 14 for the A = 3.436 double-delta planform have both been computed using 
modifications of the Multhopp approach. The curve of reference 8 is based on the two 
chordal loading function combinations used by Multhopp with a process of submatrix inver- 
sion of the matrix of influence coefficients to arrive at the local nondimensional lift and 
pitching moment. The results of CLq, and CD,i/CL^ are within 2 and 1 percent, 
respectively, of each other, whereas the aerodynamic centers differ by about 3.7 percent 
of the reference chord. The reason for this larger discrepancy in aerodynamic center 
than found for Cj_ l/y or CD,i/CL^ can be traced primarily to the different local 
pitching- moment terms computed by each procedure. The reason for the differences in 
local pitching moments is not clear. 

(5) Three span loadings are presented for the double-delta planform of figure 15 and 
are computed from the following methods: Campbell's (ref. 17), subsonic kernel function 
(ref. 18), and present. Since the subsonic kernel function and the present method are basi- 
cally identical (in the reduced frequency case) except for the handling of the integration 
across the spanwise singularity, it is not surprising that the span loadings should be quite 
similar. Consequently, the spanwise centers of pressure y cp differ only by about 2 per- 
cent semispan. Further, although the values of Cl^ are within 1 percent, reference 18 
gives a more stable pitching moment resulting in a more rearward aerodynamic center. 

Reference 17 presents a span loading which has more of its loading outboard, there- 
fore, a larger value of y cp . Also, its CLq, is lower and its aerodynamic center is 
ahead of the others. 


In figure 15(b) the same trends in the spanwise variation of the centers of pressure 
are seen to occur in each of the two lifting surface theories over the first 75 percent of 
the semispan. However, for 77 > 89 percent of the semispan, the subsonic kernel func- 
tion indicates that the chordwise center of pressure (x cp j lies ahead of the wing. The 
local center of pressure for Campbell's method is the centroid of the two-dimensional 
loading and occurs at the local c/4. 


The difficulty experienced by the subsonic kernel function approach can be attributed 

to an improper location of collocation points used in that solution. The control points 

were at — = 0.2, 0.4, 0.6, 0.8; = 0.15, 0.35, 0.55, 0.75, and the semispan of the region 

c b/2 

which contains the control point was set equal to 0.10. However, since there is not an 
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optimum location procedure specified with this method, it is unclear what number and 
location of collocation points would yield the best results. Perhaps more realistic results 
could have been achieved by locating points farther outboard than 0.75. A cosine distri- 
bution would have made an interesting comparison, but the collocation points cannot be 
located very near the tip. 

A simple case was attempted using only the outboard panel extended to the root, with 
the collocation points located the same as previously mentioned ^in terms of x/c and 

and with a region around the control point of semispan equal to 0.10. In this case 

reasonable results were found which indicate that the planform shape should be of pri- 
mary importance in selecting these points. 

Concerning the distribution of loading along the chords (fig. 15(b)), it is seen that, 
in general, both lifting surface theories have the same overall shape and differ mostly in 
magnitude until the subsonic kernel function becomes influenced by the locations of the 
collocation points, that is, when x C p is ahead of the planform. Comparing the present 
method with that of the two-dimensional distribution, the effects of induced camber can be 
seen concentrated mostly inboard of the leading- edge break. The strength of the two- 
dimensional loading of any spanwise station was chosen to be that predicted by Campbell's 
lifting- line solution. This was accomplished by requiring that the total circulation 
strength be equal to the integral of the two-dimensional distributed loading times a suita- 
ble multiplier. 

Comparisons with experimental data.- Since experimental data are generally lacking 
for the planforms previously examined, the ability of the present method to predict the 
aerodynamic characteristics of these wings cannot be appraised. Therefore, the following 
discussion will relate to three wings for which experimental data are available and com- 
parisons with theoretical results will be made. 

(1) The two predicted span loadings presented in figure 16(a) for the A = 2 delta 
planform agree closely and estimate reasonably the experimental values found in refer- 
ence 19. Also, the Cl values predicted by the two theoretical methods, present and 
vortex lattice (ref. 20), are in good agreement with the experimental result computed from 
the CLq, of the force data, since they differ by only about 5 percent. (From the variation 
of Cl with a determined from the integrated pressure measurements, the Cl value 
at a = 4.30 is lower than that found by the force tests. The difference is attributed to 
the lack of enough pressure orifices near the leading edge resulting in a poor integration.) 

A comparison of the three aerodynamic centers shows that the value given by the 
present method is ahead of both the experimental location and the one found by refer- 
ence 20. (In general, for the lower aspect ratio wings, solutions obtained using higher 
values of N and m result in a more aft location of the aerodynamic center. (See 
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figs. 6(a) and 7(a).) From an examination of the local loadings (figs. 16(b) and 16(c)), it 
is seen that although the experimental local centers of pressure are better predicted by 
reference 20, over most of the chord (x/c of 0.1 to 0.7), the present method gives better 
agreement with the experimental pressure distribution. The disagreement between the 
large pressures predicted near the leading- edge singularity and those measured may have 
resulted from the thickness effects present in the experimental data. The singularity in 
pressure at the leading edge also causes the local centers of pressure, predicted by the 
present method, to be ahead of the experimental value. 

(2) The present method predicts quite closely the experimental span loading of the 
A = 8.02 sweptback and tapered wing of figure 17 as found in reference 21. This is seen 
not only by the distribution but in the y cp location. Also, the agreement between the 
values of Cl and aerodynamic center is good. From figures 17(b) and 17(c), it can be 
seen that whereas the local chordal centers of pressure are accurately predicted, the 
pressures are, in general, underpredicted over the chord from x/c of 0.1 to 0.6 and 
overpredicted from 0.6 to 1.0. Because of the singularity at the leading edge, the theo- 
retical pressure coefficients in this vicinity are in excess of experimental values. 

(3) In figure 18, the theoretical and unpublished* experimental loading and moment 
distribution are given for a variable- sweep wing in one sweep position. The predicted 
overall Cl at ot = 3.14° is over 10 percent higher than that found experimentally from 
force data and this difference is mentioned later. The predicted aerodynamic center is 

7 percent of the reference chord ahead of the experimental value; however, this is not a 
large discrepancy when one considers that the reference chord in this case is based on 
the outer panel extended to the root. (See fig. 18.) 

The experimental values of span loading which are shown in figure 18(a) agree well 
with those predicted by theory; however, there is a slight overprediction by the theory 
over the outermost part of the outer panel which can also be seen in the c nQ! graph pre- 
sented in figure 18(b). This theoretical overprediction of local c nQ! leads to the overall 
C L being higher than the experimental value. The effect of the leading-edge shed vortex 
and different amounts of induced angle of attack generated by the theoretical and experi- 
mental loading distribution undoubtedly account for some of the disagreement seen in 
figures 18(a) and (b). Good agreement is found for the c mQ! data obtained from theory 
and experiment. The xcp graph is just the result of the Cn a and c ma curves and 
hence shows the local center of pressure slightly off of that predicted by theory. 

At the four different spanwise stations shown in figure 18(c), the chordal loadings 
are compared and the agreement is found to be reasonable at most locations. 

Illustrative ca lculations.- Effect of variable sweep: Changes in the aerodynamic 
characteristics due to increasing or decreasing the outer panel wing sweep (from the 

*Tests performed in the Langley high-speed 7- by 10-foot tunnel. 
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streamwise tip position) can also be predicted with the present method. For example, the 
planform in figure 18 has its sweep varied from 15° to 40° at a Mach number of 0.23 and 
the results, shown in figure 19, are based on the area and mean geometric chord of the 
outer panel extended to the root when A = 30°. Those results presented are from the 
present method and unpublished experimental data measured on a semispan pressure 
wing at the Langley Research Center. 

The agreement obtained between theoretical and experimental results is only fair in 
that the theory shows the proper trends and similar increments, but not absolute levels. 
Experimental values for were not available from force data. As reported 

in references 22 and 23, the variations presented in figure 19 are typical of variable- 
sweep wings. 

Effect of Mach number: The effect of increasing the subsonic Mach number on the 
wing- alone aerodynamic characteristics is accounted for by use of the Prandtl-Glauert 
compressibility rule as mentioned earlier. (See ref. 11.) In general, the results obtained 
appear to be questionable above a Mach number 0.8 because of the inapplicability of the 
Prandtl-Glauert transformation as M — 1 (ref. 11). 

Some results obtained by applying this rule to a highly sweptback and tapered wing 
are presented in figure 20 along with those determined from Mach numbers 1.2 to 2.8 by 
use of the supersonic lifting-surface theory of reference 4. Experimental values for the 
wing-body combination are presented for comparison. 

In general, the theoretical methods are able to predict reasonably well the experi- 
mental flat wing overall aerodynamic characteristics (which for subsonic speeds have not 
been published but are found in ref. 24 for supersonic speeds) even though the experimen- 
tal data contained the influence of a body. The one exception is the induced drag param- 
eter which does differ considerably. However, this is not unexpected since the real 
leading edge is sharp and the leading- edge suction is considerably less than the 100 per- 
cent assumed by the theory. Further, no induced drag, called vortex drag at supersonic 
speeds, is given supersonically since it was not available separate from the wave drag. 

In the transonic region (dashed lines) the results are obtained by fairing from sub- 
sonic to supersonic speeds. 

Warped wing: The determination of the aerodynamic characteristics for a warped 
wing can also be found by using the present method when the local surface slope distribu- 
tion is specified. In addition to the flat plate loading characteristics, which are always 
found, the basic loading features are also computed by the use of the appropriate equations 
given in appendix B. Also determined are the zero-lift angle of attack a; 0 , pitching- 
moment coefficient C m o , and root bending-moment coefficient Cg 0 . Illustrative 
results of this procedure are presented in figure 21 for a double-delta planform with 5° 
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of linear twist and 2-percent circular arc camber. (Other details of this planform and 
loading distribution are to be found in the "Supplement to NASA TN D-4427" since it is 
one of the sample cases given.) 

MEAN CAMBER SURFACE DETERMINATION 

The problem of determining the mean camber surface required to support a pre- 
scribed loading on a given planform is easier to solve from the basic formulation than 
finding the loading associated with a prescribed surface shape for the same planform. 
This can be seen from the fundamental equation in compressible subsonic flow which is 
equation (1) and is given as follows: 
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from which the mean camber surface is found by use of the following equation: 
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The determination of the local Ap(£,? 7 ) from equation (1) has been shown to involve 
solving the equations inversely for a part of its integrand; whereas, to find w/U at any 
point, the individual specified loading must be simply multiplied with the appropriate ker- 
nel function value and summed over the surface. Thus, if the more difficult problem is 
solved or programed first (as in the present study), then the other only requires, in addi- 
tion to those changes mentioned, the deletion of the matrix inversion and the addition of 
the w/U integration to find z/c. 

Since values of w/U are only computed at the control points, some representation 
at other chordwise locations is needed in order that the integration of equation (13) can 
be carried out. A least-square polynomial curve fit of the w/U values at the control 
points was chosen to represent the variation of w/U along the chord between the first 
and last control points. The polynomial curve fit is specified to have a degree of one less 
than the number of chordal control points to insure that a good fit can be obtained. Out- 
side the range of the control points, a linear extrapolation of the polynomial curve fit is 
used to the nearest chordwise edge. 
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Optimization of N and m 


As in the load determination section, a study was undertaken here to find the rela- 
tionship that N and m should have in order that the best mean camber lines could be 
obtained. It was intuitively felt that these mean camber lines would occur when the most 
control points were employed. For the A = 2 delta planform shown in figure 22, four 
combinations of N and m were selected and in each case the same surface loading 
distribution was used. It was composed of only a sin a chordwise loading which varied 
elliptically from tip to tip and was adjusted to give an overall lift coefficient of unity. 

All other pressure mode coefficients were set equal to zero. 

Figure 22(a) shows that the computed downwash w/U values at a constant m with 
varying N do not lead to a much changed polynomial curve fit, whereas at a constant 
value of N varying m does. The curves appearing in this figure are integrated and 
appear in figure 22(b). As would be expected, sets of N and m which have essentially 
the same curve fit lead to the same mean camber line. Thus, the mean camber lines for 
the sets N = 4 and N = 8 at m = ll are the same. This is also true for the sets of 
N = 4 and N = 8 at m = 23. 

Thus it appears that for the given surface loading, increasing the number of span- 
wise stations at which the desired chord loadings are defined is more important in 
obtaining the best mean camber lines than increasing the number of chordal pressure 
modes involved (much above those needed to define the shape) at each spanwise station. 
This is a similar conclusion to that reached in reference 25. (Note that the prescribed 
loading does not contain any leading- edge suction and hence it is not known what the 
effects of including the suction term would have on the camber for optimum combinations 
of N and m.) 

With the relationship that N and m should have for certain loadings established, 
some comparisons with other theoretical methods are made and follow in the next section. 

Comparisons 

The accuracy of the present method may be appraised by comparing in figure 23 
the downwash distribution at the plane of symmetry of an A = 10 rectangular wing with 
that of a two-dimensional wing with the same sin d chord loading at each span station. 

Over most of the chord the three-dimensional downwash is slightly lower than that 
of the two dimensional by a constant amount. If a much higher aspect ratio wing had been 
used in the computation, this difference would become smaller. 

At the edges, where x/c — 0 and 1, the three-dimensional downwash is seen to 
tend toward plus and minus infinity, respectively. This is caused by the logarithmic 
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singularity correction term, which results from the second-order singularity I — - — I 
in the integrand, containing a term which increases without bound at all edges. (See 7 
eq. (A9) .) 

From figure 23, it can also be noted that by using only a sin 3 loading increasing 
the number of chordal control points does not change the downwash distribution at or in 
the vicinity of the chordal control points already computed but does give a more detailed 
representation of the distribution. 

As another check, the mean camber lines obtained with the present method for the 
sweptback wing shown in figure 24 are presented there2 for comparison with those pre- 
dicted by the mean camber surface solution of Katzoff, Faison, and DuBose (ref. 26). 
Three differences are apparent in the procedures used and are 

1. The present method rounds any kinks that occur in the wing planform (i.e., at the 
root) which results in a closer representation of the isobars on the planform (ref. 5) and 
hence gives a better downwash distribution than the method used in reference 25 (ref. 27). 

2. The uniform chord loadings ACp used in reference 26 are approximated herein 
by a Fourier series of four odd sine terms - sin 3, sin 33, sin 53, and sin 73. This 

Q 

series or others like it composed of cot ^ and even sine functions give greater flexibil- 
ity in approximating arbitrary pressure loadings. The amplitudes of the chord loading 
were uniform over the entire area and produced a lift coefficient of 1.0. 

3. The present method uses the least- square polynomial curve fit and constant slope 
scheme to determine the downwash values all along the chord as discussed earlier. This 
procedure differs from that used in reference 26 where the downwash distribution com- 
puted at selected internal constant chord lines is completed by fairing except near the 
edges of the chord. There the downwash is assumed to be equal to that required to pro- 
duce a constant pressure all along the resulting mean camber lines. 

The curves which are discussed in the previous paragraph are then integrated to 
find the mean camber height z/c above the x-y plane. Despite the differences dis- 
cussed between these two methods, figure 24 shows that the agreement between them is 
good except near the plane of symmetry. This is explained in part by the first differ- 
ence listed. (See also ref. 27.) 

Application .- An application of the present method to the highly sweptback and 
tapered planform shown in figure 25^ is made where the lift coefficient for this planform 

^The results of reference 26 were found by integrating from the leading edge rear- 
ward and hence were adjusted in figure 24 to the same zero level for comparison. 

^Tabular results for this planform and loading distribution are to be found in the 
"Supplement to NASA TN D-4427" since it is used as the sample case. 
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was taken as unity at a Mach number of 0.30. The chordal loading function coefficients 
were specified so that there was (1) a zero cot distribution, (2) an elliptical distribu- 
tion of sin 3 across the span, (3) a sin 23 distribution such that the local x c p varied 
parabolically from 50 percent of the local chord at the plane of symmetry to 20 percent at 
the break and from the break linearly back to 50 percent at the tip, and (4) a sin 33 dis- 
tribution which caused the slope of the chordal loadings to be zero at the trailing edge. 

Figure 25 shows that the mean camber lines have a rearward movement of the x/c 
location of maximum camber and a twist which increases toward the tip. The mean cam- 
ber line at ?7 = 0.997 is typical of stations near the tip and should be examined carefully 
before a direct application is made since the downwash distribution which produced it may 
be overly influenced by the singular nature of the downwash at the tip. 

CONCLUSIONS 

The results obtained with a modified version of Multhopp's lifting- surface theory 
programed for the IBM 7094 electronic data processing system and Control Data 6400 
computer system to find the surface loadings and mean camber surfaces have been 
studied and the following conclusions are drawn: 

For the surface loading program: 

1. It has been found that the loading program may be applied with reasonable accu- 
racy to conventional and composite planforms to find their flat plate characteristics and to 
wings with twist and/or camber to determine the basic loading distribution as well as the 
force and moment coefficients at zero lift. 

2. In predicting the surface loading, the number and spacing of control points is more 
important when seeking an accurate overall aerodynamic center than when a solution for 
the lift- curve slope, span loading, or spanwise location of center of pressure is being 
sought. It has been found that at least four, and preferably more, chordal control points 
should be used in obtaining aerodynamic centers. This number should depend upon the 

/3A of the planform being considered, since the more control points used (100 available) 
at each spanwise station necessitates a decrease in the number of spanwise stations. This 
may lead to a problem when values of j3A are large. However, in general, the number 
of chordal control points should be kept between 6 and 10 to insure that the resulting val- 
ues of aerodynamic center lie within the region of converging results. 

3. Reasonable agreement is generally obtained with both experiment and other theo- 
ries concerning local and overall loading characteristics. 

4. The effect of Mach number can be determined and the resulting agreement with 
experiment is found to be reasonable. 


For the mean camber surface program: 

1. Good agreement is achieved between the downwash required by two-dimensional 
theory for a sin 3 loading and that predicted by a three-dimensional uniform sin 3 
loading at the plane of symmetry of an aspect-ratio- 10 rectangular wing. 

2. The agreement with results from other mean camber predictions is found to be 
fairly good and the present method has the capability of approximating most chord loading 
shapes by the use of a Fourier sine series. 

3. For an acceptable prediction of mean camber surfaces which includes no leading- 
edge suction, the number of chordwise stations should be nearly equal and never less than 
the number of pressure loading terms required to properly define the desired shape and 
the number of spanwise stations should be as many as can be accommodated by the 
machine for that number of chord stations. 

4. The resultant mean camber surface predicted for a composite planform whose 
loading includes no leading- edge suction was found to have a reasonable twist and camber 
distribution inboard of the tip. In the vicinity of the tip, the mean camber surface does 
not appear to be as reasonable due to the singular nature of the downwash at the tip. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 19, 1968, 

126-13-01-50-23. 
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APPENDIX A 


FORMULATING THE INFLUENCE COEFFICIENT MATRIX 


The formulation of the influence coefficients L^ n (0s) is accomplished by integra- 
ting (or approximating) the chordal influence that the pressure loading functions exert on 
each control point from every other spanwise station. (See ref. 6.) Then, for symmetri- 
cal spanwise loadings, the chordal influence at r]_ n> L],_ n (0 s ), is added to that of ri n , 
Lj, n (0s)> to give the total L^^g), and for antisymmetrical loadings the influence at 
rj_ n , L|,_ n (0 s ), is subtracted from that at rj n , Lj m (^ s ). The value at the plane of sym- 
metry, L^s), is just equal to l|o(^ s ) when the loading is symmetrical and zero 
when antisymmetrical. 


Now, for v * n, 


with 


v 


kin($s) = -bz/nKj($s,y ulVn) 
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and the subsonic kernel function is given by 

£c(y)cos <fi + d(y) + c{t])cos 3 - 


K(<p,y;3,r)) = 1 + 


^ c(y)cos (p + d(y) + c(?])cos 3 - dO?)] 2 /# 2 + (y - n /) 2 


(A4) 
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where 


x = -c(y)cos <p + d(y) 
£ = -c(?7)cos $ + d (77) 


(A5) 


Now, for n = v, 


where 


Kj(4> s ,y^) = \ 


+ Logarithmic singularity correction 

m + 1 

vv " 4 sin 9 V 


r 

2(p s + 2 sin 0 S 

0 s 

II 

*33 

(p s ' 2 Sin 2( f )s 

(j = 1) 

sin(j - l)<p s sin(j + l)<j> B 

(j >1) 

j - 1 j + 1 


(A6) 


(A7) 


> (A8) 






The logarithmic singularity correction term (LSC) arises because of the spanwise inte- 

■ 1 

Its derivation is discussed by 


gration across the second-order singularity . 

(y - vr 

Multhopp in reference 5, Mangier and Spencer in reference 28, and Van Spiegel and 
Wouters in reference 6. Its best form is 


LSC = 




(dfj(0)/d0) 0=0s 

m + 1 

c(y y ) 2 (m + l)sin Q v sin <£ s 

8 




-(cos 29 v - In 4) 


m-1 

2 




sin2fl n In |yj, - r? n | 


(In - v\ odd) (A9) 


where the prime on the summation sign means that the term when n = v is eliminated in 
the summation process and 
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*iW 




(A10) 
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EQUATIONS USED IN THE COMPUTATIONS 


The following equations are used to determine the total local characteristics at the 
design lift coefficient as well as the overall characteristics for wings with twist and/or 
camber. Note that the qtj(? 7 n ) terms represent the pressure mode coefficients asso- 
ciated with the twist and/or camber boundary conditions, and the qj(?7 n ) terms represent 
the pressure mode coefficients associated with the flat plate boundary conditions. The 
qtj ( Vn ) terms are arrived at in a manner similar to that of the qj (?7n) terms, as given 

-4w(0 s ,y y ) 

by the matrix inversion of equation (4), except that the ^ values are now the 

local angle of attack for the warped wing at the flat plate control points when the root 
chord is at zero angle of attack. 

Local Aerodynamic Characteristics 

(a) Pressure coefficient 


AC p O,„ n ) = -j-i-iLhn) - 

p c (%)qoo ] \ c L,q!= 1 rad / 


cot i- 


^L,l “ C L, desired] 


+ l_ n (%) • q i (,?n) r"cz^r 

i=i 


rad 


sin j3> 


j 


(b) Chord loading 


( c l cq oo) n = W.qtofon) + \ “ l^o^n) 


1 . ^|/Cl,1-C L , 

c L , a=1 



(c) Nondimensional lift coefficient (called circulation in the "Supplement to NASA 
TN D-4427") 


y n - ( 2 b ) n ‘ 2^ oo ^ t O (77n) + 2 qtl(77n ^ 
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(d) Pitching moment about local leading edge 


(c m c2 qoo ) n = c(? 7n ) 7 T^t 0 (r ?n ) + - |qt 0 (? 7 n ) " \ 

-EoW - ! <!,(»„> - §q 0 W - Zhems* 

*— * J \ '-'L,q!= 1 rad 


(e) Location of center-of-pressure loading from Y-axis 

,Y C L,1 " desired' 


(^cp) n “ ” 


c ( Vn ) j qto^n) + | - q 0 (? 7n) + | q2^n) 


Cl,q!= 1 rad 


fetolvn) + iqtl(Tta)] - [q 0 <»n> 


) + d (Vn) 


(f) Center-of-pressure loading in fractions of local chord 

(?cp) n “ d (v n ) + C <V 


( x cp) n = 


2c(rj n ) 


(g) Span load coefficient 


= * 


\ Ca % c ayqoo 


tfofan) + 1*1*1 fan)] ' [^ 0 (r ?n) + \ ^l^nj] ^ 


desired\ 
rad 


Overall Aerodynamic Characteristics 

(a) Flat plate lift 

m-1 

2 

I J 1 H^oo dr l * y (50^11) + \ qi(^? n )] sin e n 

n=-S^i 
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(b) Flat plate pitching moment about the total wing c/4 


m-1 

2 


!l 1 1 cmc2q ” 


dr/ 


7r2(b/2) 
m + 1 


k) - x c - fJj^O^n) + 


n= 


m-1 


c(?7n) 


| + \ <l2^n)J> sin °n 


c ref 

Xc ref " “4 


(c) Flat plate pitching moment about the reference wing Cref/4 = 

(Flat plate pitching moment) total wing - Lift^cg- + J 

(d) Flat plate rolling moment about the X-axis 

m-1 

2 

(I ) 2 I 1 V ° l cq “ 1 / + | qi<?Jn)] sin ®n 

-v 

(e) Flat plate location of the aerodynamic center 


( ?c p), 


(Flat plate pitching moment) reference wing c ref 


overall 


Lift 


+ + 


Cref 4 


(f) Flat plate aerodynamic center in fractions of reference chord, measured from 


the x 


c ref 


a.c. = 


(Flat plate pitching moment) re f erence w i n g ^ 


Lift(c r ef) 


+ 4 


(g) Angle of attack at zero lift 


oi n = 


C L,1 

C La 


where Cl ,1 is determined similarly to flat plate Cl 
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(h) Pitching-moment coefficient at zero lift 

Cm,o = Cm,l + ^m a a o 
where C m> i is determined similarly to flat plate C m 

(i) Root bending- moment coefficient at zero lift 
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(a) Planform variables. 



(b) Reference wing geometry. 

Figure 1.- General planform which can be used in these programs along with its coordinate system and reference wing geometry. 







Number of span stations, m 

(a) Induced drag parameter. 

Figure 4.- Effect of varying number and locations of control points and chordal loading functions on aerodynamic characteristics of A = 2 rectangular planform at M = 0. 









(a) Induced drag parameter. 

Figure 5.- Effect of varying number and locations of control points and chordal loading functions on aerodynamic characteristics of A - 7 rectangular planform at M = 0. 












Number of span stations, m 


(a) Aerodynamic center. 

Figure 6.- Effect of varying number and locations of control points and chordal loading functions on aerodynamic characteristics of A = 1.7 double-delta planform at M = 0 







Number of span stations, m 


la) Aerodynamic center. 

, double-delta planform at m - U- 

and locations of control points and chordal i.ding functions on aerodynamic character, ,cs 


Fioure 7- E«ect of varying number 









Figure 8.- Span loading coefficient variation with different combinations of N and m for A = 2 rectangular pianform at M 


E 







« >oo □ 



(c) Constant m. 
Figure 8.- Concluded. 
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Figure 10.- Span loading coefficient variation with different combinations of N and m for A = 1.7 double-delta planform at M - 0. 












I 



Number of span stations , m 

(a) Aerodynamic center. 


Figure 12.- Effect of varying number and locations of control points on aerodynamic characteristics for four A = 2 planforms at 

M = 0 when N = 4. 
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(b) Local spanwise centers of pressure and chordwise load distributions. 
Figure 17.- Continued. 



N= 4 , m=37 



I .2 .3 4 .5 .6 7 

x/c 

(b) Concluded. 


Figure 17.- Concluded. 




(a) Spanwise. 

Figure 18.- Loading and moment distributions predicted by present method and compared with experiment for A = 4.303 variable- sweep 

planform at M = 0.23. 









Figure 20.- Effect of Mach number on some longitudinal 




Figure 21.- Warped A = 1.491 double-delta planform at M = 0 and computer results. 
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A -10.000 




0 ./ 


— Two-dimensional downwash due to uniform 

sin 0 chord loading 

— Three-dimensional downwash due to uniform 

sin 0 chord loading ( A - 10 rectangular wing} 


Figure 23.- Comparison of two- and three-dimensional downwash predictions from uniform sin 9 chord loading at n = 0.000. 






National Aeronautics and Space Administration 
Washington, D. C. 20546 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS A! 
SPACE ADMINISTRATION 


OFFICIAL BUSINESS 


FIRST CLASS MAIL 


0 10 OUl 2 - 3 l)S 6b 10^ 00903 

A l % i- O - v C L 'a fc A u N S LABL^ATu'O/^f h L / 

K I < i L A •; 0 L K P i .IK C l B A S l- t * 1 ^ *" 1 ^ I 0 u b / I 19 


A I 1 M I N S '■> U F L I K f- F - v, ^ \0 \t A y 

L 1 H A ?< V /v L IL / 


u • * 1 1 J l t L i : I 


POSTMASTER: 


If Undeliverable (Section 15 
Postal Manual) Do Not Ren 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of inf or ination concerning its activities and the results thereof .” 

— National Aeronautics and Space Act of 1958 


NASA S CIEN TIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. • 

TECHNICAL NOTES: Information less broad 
in scope but 'Nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C 20546 


